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ABSTRACT 


A meriu driven int^er active Pascal program has been developed 
for the thermal hydraulic design of shell and tube heat 
exchangers. The manufacturing tolerances and data for 

geometrical clearances are taken from THEMA for class R 
exchangers. The program can be used for rating as well as for 
the design purposes. The Kern method is used to get the primary 
design which is further modified by Bell method. Finally the 
checks of possible vibration damage because of vortex shedding 
and fluid elastic whirling in addition to the baffle and 
collision damage numbers are carried out. 

The results are encouraging when checked for the avilable 
case studies* considering the complexity of the problem. 


NOMENCLATURE 


Subscripts s»t and m are used for shell side* tube side and 
tube metal properties respectively. 

A Surface area> 

2 

Ab-a Bypass area of one baffle*. m 

2 

Amh Minimum flow area at center line of one baffle*. m 

2 

A«b Shell to baffle leakage area* m 

2 

Atb Tube to baffle leakage area* m 

2 

Aw Wet flow area through one window* m 

2 

Awg Gross window flow area* m 

2 

Awt window area oooupied by -tubes, m 

Be Baffle cut, ?S of shell diameter 

Bt Baffle thickness, m 

Cn Frequency constant 

Cp Specific heat for constant pressure, J,-^Ckg.KI> 

Cv Specific heat for constant volume, J/^Ckg.O 

Dctl Pitch Diameter of outer most tube, mm 

Db Baffle Diameter mm 

IX Inside diameter of tube, mm 

Doll Outer tube limit diameter, mm 

E>s Shel 1 Di ameter mm 

Dt Tube outside diameter, mm 

Dw Equivalent Hydraulic diameter of baffle window, mm 

Modulus of elasticity of tube metal , 


E 



Fc 

Fv 

f 

fr, 

f^va 

S 

Gw 

gc 

hv 

ha 

I 

j 

Jc 

Jl 

Jb 

Jr 

Js 

k 

L 

Lb 

Lbb 

Lp 

Lpn t Lpp 

Lcb 

Ltb 

Lts 

Nb 


Fraction of tubes in cross flow 

Fraction of tubes in one window 

Friction factor 

Nat ur al f r equency » Hz 

Vortex shedding fr equency » Hz 

S£ 

ma.ss v^&l oei -t y ^ k g^'C m see . _> 

Wi ndow mass vel oei "ty kg/C m see . I? 

Gravi National conversion eonslant* Cin S. I. unit, = 13 
Ideal Bundle hea*t ’transf’er eoef’f'i ei ent, CW^‘m*lC3 
Cor recited shell side heat^ t,ransf'er eoef^fiei ent, 

Moment, of i ner t, i a » 

Heat transfer f actor > given in appendix 1 
Correction factor for baffle cut» to calculate ha 
Correction factor for baffle leakage effects 
Correction factor for bundle by pass effects 
Correction factor for adverse temperature gradients 
Correction factor for variable baffle spacing 
Ther mal conduc t i vi ty » C W/mK3 
Length * m 

Baffle spacing C Lbi = inlet* Lbo = outlet* Lbc = 
central 3* mm 

Bundle bypass diametral gap* mm 
Bypass lane* mm 

Tube pitch dimensions in appendix 1 
Diametral clearance* shell to baffle mm 
Tube to baffle diametral clearance* mm 
Thickness of Tubesheet* 

Number of baffles 


mm 





(Mi> 

Number of pairs of sealing strips 
Nt Total number of tubes in exchanger 

Ntcc Number of effective tube rows in cross flow 

Ntcw Number of effective tube rows in baffle window 

Rh Correction factor for bypass flows for for AP« 

eaicul ation 

Ref Reynold Number C subscript « = shell > t = tube sideZ? 

Rt Correction factor for leakage effects for AP« 

calcul ation 

Rs Correction factor for inlet and outlet zone flow 

expansion 

Sic Strouhal number » 

T Temper at ure> "^C 

U Overall heat transfer coefficient^ W^'Cm KZ? 

Vc Cross flow velocity* nru^sec. 

VcrLt Critical cross flow velocity* m/sec. 

Greek syirxibols s 

ft Instability constant used in Conors equation 

Logarithmic decrement of amplitudes because of damping 
Dynamic viscosity of fluids* N. s/m^ 

Dynamic viscosity of fluid at wall temperature 
p Etensity kgx"m^ 

& Angles defined in fig. 5* radian 



CHAPTER 1 


INTRODUCTION 


The use of computer programs is essential for thermal 
hydraulic design of heat exchangers as they enable to investigate 
rapidly all parameters to arrive at an optimum solution. 
Manufacturers use computer programs to design heat exchangers 
that will meet the required duty within certain operating 
constraints. Users, on the other hand use programs to check what 
a given exchanger will achieve if operated in a specified way. 

This thesis is concerned with the development of a computer 


code for thermal hydraulic design 

of 

shell 

and 

tube 

h&a'L 

exchangers. A brief introduction 

of 

shell 

and 

tube 

heatr 


exchangers including their types, and criterion for selection and 
design are given in following sections 1.1 to 1.3. 

1 . 1 SIhell and Tube Heat Exchangers, Ck^iqponents and Designation : 

Heat Exchangers built of round tubes mounted in cylindrical 
shell with their axes parallel to that of the shell are called 
Shell and Tube Heat Exchangers. It is the most common of the 
various types since it accounts for over 50>S of all installed. 
Although it is not compact, it is robust and its shape makes it 
well suited to pressure operation. It is also versatile and it 
can be designed to suit any application. Except for the special 
purpose air cooled heat exchanger, it is usually the only type 




1. stationary Head — Channel 

2. Stationary Head — Bonnet 

3. Stationary Head Flange — Channel or Bonnet 

4. Channel Cover 

5. Stationary Head Nozzle 

6. Stationary Tubesheet 

7. Tubes 

8. Shell 

9. Shell Cover 

10. Shell Flange — Stationary Head End 

11. Shell Flange — Rear Head End 

2. Shell Nozzle 

3. Shell Cover Flange 

4. Expansion Joint 

5. Floating Tubesheet 

6. Floating Head Cover 

7. Floating Head Flange 

8. Floating Head Backing Device 

9. Split Shear Ring 


20. Slip-on Backing Flange 

21. Floating Head Cover — External 

22. Floating Tubesheet Skirt 

23. Packing Box 

24. Packing 

25. Packing Gland 

26. Lantern Ring 

27. Tierods and Spacers 

28. Transverse Baffles or Support Plates 

29. Impingement Plate 

30. Longitudinal Baffle 

31. Pass Partition 

32. Vent Connection 

33. Drain Connection 

34. Instrument Connection 

35. Support Saddle 

36. Lifting Lug 

37. Support Bracket 

38. Weir 

39. Liquid Level Connection 


\ 



Flg.1. HEAT EXCHANGER, CONSTRUCTION 
AND TERMINOLOGY (TEMA). 
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which can be considered f‘or large surface areas having pressures 
'than 30 Bars and tamper atur-&s greater than 360 d. 

The nomenclature typically used in practice and Recommended 
by Tubular Exchanger Manufacturers Association^ C TEMA D for heat 
exchanger components is shown in Fig.l. 

Itony variations of shell and tube exchanger type are 
aval 1 abl e > depend! ng upon the combi nati ons of shel 1 i and head C 
front and rear end head 5 types. The difference lies mainly in 
the detailed features of construction and provisions for 
differential thermal expansion between tubes and shell. The 
Shell and head types are described in section 1.1 to 1 . S. 

1.1.1 Head Types 

TEMA classification for rear and front end heads is as shown 
in Fig. 2. and described as under. 

1.1. 1.1 Channel type head CTEMA A,L.C»N5 

This type of head is used when cleaning of inside of tubes 
is to be performed frequently, because it permits access to tube 
si de wi thout di stur bi ng pi pi ng connect! ons . It can be ei ther 
integral with shel 1C TEMA C,N3 or bolted to it CTEMA A.L>. 

1.1. 1.2 Bonnet type head CTEMA B,M:> 

It is usually fitted at rear end of fixed tube sheet 
exchangers when nozzles are not there and cleaning of tube side 
is not frequent, because it permits access to tubes only after 
removing piping connections, but it is cheapest design. 

1.1.2 iShel 1 Types 

Depending upon nozzle arrangement on shell side, there are 
seven TEMA standard arrangements as shown in the fig. 2. ITrie most 



5HEU TYPIS 


HEAR END 
HEAD TYPES 




Fig.2. SHELL, FRONT END. AND REAR 
HEADS. TEMA CLASSIFICATION 
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common arrangomen-t is E t<ypa sholl which has one shell side pass. 
Other arrangements are F type which has two shell side passes and 
inlet and outlet nozzles on same end of shell is used^ when it is 
needed to increase shell side flow velocity hence pressure drop 
with out increasing the length of shell. The split flow 
arrangements Ctype G, J!> are used to reduce shell side 

pressure drop. Type K is used in kettle type reboilers used to 
vaporize the shell side fluid, and cross flow shell type X also 
has low pressure drop characteristics Clower than J type!) because 
there are no segmental baffles so shell side fluid makes only one 
pass across the bundle. 

In horizontal units inlet should be at the bottom and outlet 
at the top rather than both at the sides. This minimizes the 
bypass of shell side fluid. 

1.1.3 Desi gnatl on 

TEMA Recommends that Heat exchanger typie and size be 
designated by shell inside diameter in inches rounded to nearest 
integer followed by tube length in inches. For example, 

SIZE 19-S4 TYPE BGU is the designation for a U tube 
exchanger C TEMA ID with bonnet type stationary headCTEMA BD , split 
flow shellCTEMA G5 , 19" inside diameter shell with tubes 7* 

straight length. Classification written in parentheses is that 
described by TEMA ,and is shown in Fig. 2. 


B 
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1.2 Principal -types <at shell and ttdse heat exchangers 

1.2.1 Fixed Tube sheet CTEMA L, N:> 

As the name Implies tube sheets are welded to shell > so 
thermal expansion is the criterion f^or its selection. Although 
it is an extremely leak proof design but shell side cleaning can 
only be done by chemical washing. Expansion Joints in shell, if 
provided, reduces over -stressing of components because of 
differential thermal expansion. It is the least cost and easiest 
to fabricate type. 

1.2.2 U tube CTEMA UD type 

It needs only one tube sheet, and U tubes are free to expand 
relative to shell so problem of thermal expansion is eliminated 
and design is suitable for high pressure and temperature. Tube 
bundle is removable so mechanical means can be adopted to clean 
on shell side, but tube side cleaning becomes difficult because 
of bend. For clean tube side fluids with higher pressure and 
temperature opjorations this is ideally suitable. 

1.2.3 Floating head type 

In this type one tube sheet is fixed but other is free to 
float. There are 4 different types under this category, which 
are described below. 

1.2. 3.1 Split backing ring CTEMA SD 

To separate the shell and tube side fluids at floating head 
end, the tube side of the floating tube sheet is fitted with a 
flanged, gasketed cover at its periphery which is held in 
position by bolting it to a split backing ring on the other side 
of the tube sheet. Tbe complete floating head assembly is 



7 


locaL-ted beyond bhe main shell in a shell cover of larger 
dlame-ter. The nabure of floating head eonsbruction lind'ts it, t,o 
an internal Ctube sidel> design pressure of the order of 7SO 
Ib^^in* CSO bar!>. Outer tube limit C OTL is limited by the 
inside diameter of floating head gasket or inside diameter of 
backing ring whichever is smaller. 

1.2. 3. 3 Pull through CTEMA T) 

In this type floating head cover is directly bolted to 
floating tube sheet so its diameter is larger* which further 
increases the shell diameter approximately the same as the 
enlarged shell cover of the split backing ring type to contain 
same number of tubes. The important advantage is easier 

dismantling and tube bundle removal reducing downtime and 
maintenance. Shell bundle clearance is largest of all types. 

1.2. 3. 3 Packed latern ring type or externally sealed type CTEMA 
W> 

The separation of shell and tube side fluids at floating 
head is obtained by means of packing rings installed between the 
outside of the floating tube sheet and recesses in the rear head 
flanges. The shell and tube side fluids each have their own 
packing rings* which are separated by a latern ring provided with 
weep holes for leak detection. It is the least expensive 
removable bundle type exchanger. Generally used below 200 *C and 
300 psi . Only one or two tube passes can be used. 

1.2. 3. 4 Outside packed CTEMA P:> 

In this type shell side fluid is still contained by glanded 
packing as in externally sealed CTEMA VD type but tube side fluid 



is ee]rnt,alned by metallic seal of machined surfaces so there is no 
restriction of pressure^ temperature and corrosiveness of tube 
side fluid. 

1.3 Thermal l>eslgn Parameters and Selection Criterion 
1.3.1 Tube Di ameter 

The selection of tube diameter depends upon fouling nature 
of fluid* space available and cost. Compact and economical units 
however are obtained by using small diameter tubes. Tubes of 3/4 
and 1 in. are most widely used. 

1.3.3 Tube Thickness 

The tube thickness is designated by Bermingham Wire Gauge C 
BWG Z>. The standard thicknesses for different nominal diameters 
of tubes are given by TEMA. The minimum tube thickness is 
selected considering the pressure differential across the wall* 
corrosiveness of fluids, flow Induced vibration. Axial strength 
and standardization. 

1 . 3. 3 Tube Length 

Tubes are selected as long as pxissible considering 
constraints of handling and space available at the time of bundle 
'emoval , because cheapest construction with given surface area 
dll have smallest shell diameter. The maximum tube length for 
emovable bundle exchangers is about Q m with bundle weight 
.pprox. 20 tones. For fixed tube sheet exchangers the tube 
ength can be up to IS m. Tube lengths of 3,12,16,20 and 24 ft 
re used as standard lengths and recommended by TEMA. 



a 



Fig. 3. TIE RODS AND SPACERS 


lO 

1.3.4 Tube Pitch and Pattern 

Use of small pitch Ctube to tube center distanced? results in 
compact exchangers j but again this is limited by fouling nature 
of shell side fluid and tube to tube sheet welding technique 
used. 

Square and diamond patterns of tubes permits external tube 
surface cleaning by mechanical tools* but triangular patterns 
will not so use is limited to clean service. Standard pitch to 
tube diameter ratios in use are 1.36* 1.33 and 1.5. Triangular 

pattern results in larger pressure drop and heat transfer 
coefficient on shell side. 

1 . 3. 5 Number of Tube Passes 

Multiple tube passes are used to increase the tube side 
velocity so as to lower the rate of fouling and increase heat 
transfer coefficient within allowable pressure drop. Usually 
even number of tube passes are used because of mechanical 
difficulties associated with odd tube passes in providing tube 
fluid headers on both sides of exchangersC particularly floating 
head typesD. Generally up to & tube passes are used* because to 
provide passes the head is fitted with fiat metal plates* C known 
as partition plates* shown in Fig.3. which are welded to the 
tube sheet. This makes design less compact and difficult to 
access. Tube passes are arranged in such a way that the tube 
side can be drained and vented. 

1.3.6 Sped al Tubes 


Bi-metallic tubes are used when single metal is not 
compatible with both shell and tube side fluids. Low fin tubes C 
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in which t,he fins af nominal height, l.S mm are extruded from the 
parent tube D are considered if shell side heat transfer 
coefficient is to be increased. Sometimes tube inserts are used 
if tube side heat transfer coefficient is to be increased. 

1.3.7 Baffles 

The use of baffles in shell and tube heat exchangers 
increases cross flow velocity thus heat transfer coefficient and 
reduces fouling. Baffles also act as tube supports against 
sagging and vibration. Different types are as follows- 

edifice Baffles are circular disks having holes larger 
then tube outer diameter which permits the shell fluid to pass 
through. These can be used only with very clean shell fluids. 

CS'J Disc and E)oughnut are essentially like orifice baffles 
but are cut into two piortions a disk and other is called 
doughnut . 

Segmental Baffles are having baffle holes smaller than 
in orifice and disc - doughnut types. 

The clearance between baffle hole and tube outer diameter is 
Kept as small as possible so that tubes can be made to pass 
through. Segmental baffles can be single segmental > double 
segmental or triple segmental. Single segmental baffles give the 
ttaximum cross flow velocity»heat transfer coefficient and 
pressure drop. Double and triple segmental baffles are used when 
Hsrmissible pressure drop on shell side is smaller and designer 
an not afford to have single segmental baffles. 

The cut portion of single segmental baffle is called window 
nd sp>ecified by of shell inner diameter. Smaller the baffle 
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cut, more will be the amount of cross flow but pressure drop will 
also be more, also larger baffle cut will result in 
maldistribution of shell flow and small heat transfer 
coefficient. Baffle cut of 20 to 2S is considered optimum. 

Generally baffles are spaced at a distance smaller than or 
equivalent to shell diameter and more than one fifth of the shell 
diameter Cref.TEMA3. Baffle spacing influences the shell side 
flow velocity so lower limit for baffle spacing is limited by 
ma>d.mum allowable shell side pressure drop and performance loss 
resulting from increase in leakage streams. Upper limit for 
baffle spacing is due to reduction in heat transfer coefficient 
because of increase in parallel flow and increased susceptibility 
to flow induced vibration because of decreased stiffness of tubes 
hence less natural frequency. Distance between the adjacent 
baffles is generally kept equal except at the ends where it is 
generally more than central baffle spacing to accommodate inlet 
and outlet nozzles and impingement plate etc. 

1.3. B Routing of Fluids 

Allocation of fluids to shell or tube side dep>ends UF»on 
several factors £.e. 

CID Corrosive nature of fluid : - 

more corrosive of the fluids is to be placed on tube side 
because it is easier to replace damaged tubes individually. 

C2I> Fouling nature : 

Shell side is difficult to clean than tube side except in 
he U tube exchangers, so fouling fluid is routed to tube side. 
C33 Fluid Temp»erature and Pressure : 
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It- is economl cal and safer i-o eont-ain t-he high pressure and 
t-emperat-ure fluid in t-ubes t-han in shell. 

C45 Allowable Pressure l3rop : 

For same amount- of pressure drop higher heat, t-ransfer 
coefficient, is obt-ained on t-ube side^ so fluid wit,h less 
allowable pressure drop is placed t-o t,ube side. Also prediction 
of pressure drop on the tube side is more accurate, so when it is 
crucial to limit pressure drop for a particular fluid, it is 
placed to tube side. 

CS5 Flow rate : 

Because Reynold* s number for turbulent flow on shell side 
is around 200, so smaller flow rate on shell side will experience 
more turbulence. On the other hand larger flow rate on shell 
side will require multiple tube pass construction with consequent 
temperature drop efficiency, so better overall design is obtained 
when smaller quantity fluid is placed on shell side. 

C7D Hazardous or expensive fluid is to be placed on tightest 
side which is generally tube side. 

CS!> In the case of no restrictions as listed above. The 
arrangement which gives higher overall heat transfer coefficient 
with out exceeding allowable pressure drop is chosen. 

objectives w study 

The main objective is to develop a simpler code than 
proprietary codes for the thermal -hydraulic design of shell and 
.,ube heat exchangers, which can also be used for the rating and 
"low induced vibration checking of existing heat exchangers 
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wi-Lhout, the requirement of" lerge memory and computing system. 
The software developed under this thesis is best suited to be 
used with Personal Computers » with at least 256 K RAM. 



CHAPTER 2 


DESIGN METHODS AND LITERATURE SURVEY 


2.1 General Classlflca-tlon 

Thr-ae different, t.ypes of approaches discussed in lit,erat,ure 
are as follows 

2.1.1 Int,egral Approach 

Ttiis approach considers total flow as effective in heat 
transfer and pressure drop calculations and introduces simple 
correction factors. The procedures based upon this approach are- 

C A St 3 C 

Kern*a method » and EJonohue's method 

2.1.2 Semi Analytical Approach 

This approach takes in to account the affect of non 
idealities and uses correction factors based upon the estimation 

C dL 3 

ot diitt&r&nti Gt^r^a.nAS. m^*thod is of 

this approach which is discussed in de'tail la*ter. 

Sal. 3 St^ream Analysis Metxhod 

C SE<43 

This is the most complex approach proposed by Tinker. 

[tiis method is the basis of Heat Transfer Research Institute C 
ITRI ]> and Proprietary Computer programs. Several researchers 
tave proposed their simplified methods based upfon this approach. 
Since later discussions refer to the various flow streams 

I 24' 3 

s dottignat^od by Tinkor , a brief oxplana*tion is as follows- 
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Tink&r divided t,he shell flow int/O four strr earns. As shown in 
fig. 4A. A is -the leakage st.ream -through the tube-baffle 

clearance in one baffle, B is the cross flow stream through the 
bundle, C is the bypass stream between the shell and the outside 
of the tube bundle, and E is the leakage stream through the 
baffle-shell clearance. An F stream C not shown in fig. 5 is also 
defined for multiple tube pass exchangers, because of flow along 
the pass partition plate. 

The relative effectiveness of the flow streams towards heat 
transfer is B> A> CG, F3>E. so C, F, E and A streams should 
be minimized in comparison of B stream. It is accomplished by - 
C1Z> Using small manufacturing clearances. 

C2D Using Sealing strips to minimize C stream particularly 
in floating head designs. 

C3D By avoiding too small C < 0. 2»Shell ID D baffle spacing. 
C4:3 By avoiding laminar flow on shell side, C,F,and E 
streams can be minimized. 

CSD Fraction of B stream reduces with reduction in shell 
diameter, so careful designing is needed. 

3.2 Bell method 

The following section formulates the method, all equations 
ire given in consistent units so no conversion factor is needed 
'or S I units. The nomenclature is also given but unit is given 
it the beginning with complete nomenclature of all symbols used. 

A simpler method was developed by Bell under the Delaware 
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Universi-ty Research Program' and furt-her refined by him and 

£ J 

o-t her modif'iea-tiorxs and impr-ov€?men'ts hav<& been naade by Tabor ak 
mainly in the Rh and Ri curves and some changes in Jb and Jl 
daba. 

The method provide calculation procedure Tor shell side heat 
transfer and pressure drop. 

The approach first calculates the ideal heat transfer and 
pressure drop on the basis of total flow across the bundle and 
then five correction factors are applied to approximate the 
effect of leakages and bypassing. 

The basic equation is 

h« = hi ^ Je. J t. Jb. Jtt. Jr CS. 1I> 

where hi = Heat transfer coefficient for pure cross flow in an 
ideal tube bank assuming the entire shell side stream 
flows across the tube array. 

J<s = correction factor for baffle cut. This correction 
includes the effect of window and bundle heat transfer. 
It has the value 1.0 for an exchanger with no tubes in 
window and increases to l.iS for small baffle cuts and 
decreases to 0. for large baffle cuts. For a typical 
well designed heat exchanger the value is near 1.0. 

Jl = correction factor for baffle leakage effects including 
both tube— baffle and baffle-shell leakageCA and E 
streams!? with a heavy weight is given to latter. It is 
a function of the ratio of leakage to cross flow area 
and a function of the clearances. A typical value of Jl 
is in the range of 0.7 to O, S. 


IQ 

Jb = correcbiori factor for bundle bypass CC and F sbreamaO. 
For a fixed tube sheet construction Jb O. Q and for a 
puli through Jb O. 7 which can be increased to O. Q by 
the use of sealing strips. 

J« = correction factor for variable baffle spacing at the 
inlet and outlet sections. The nozzle locations result 
in larger spacing and lower velocities and thus lower 
heat transfer coefficients. Ja usually ranges from 0. 8S 
to 1 . O. 

Jr = correction factor for any adverse temperature gradient 
buildup in laminar flows. This correction applies only 
for shell side Reynolds Wumbers of less than lOO. 

Shell side pressure drop is the summation of pressure drops 
in the inlet and outlet sections and the baffle central section* 
consisting of the cross flow and window pressure drop. The 
following correction factors are used. 

Rl = Correction for leakage C A and E streams!). This correction 
is based on the same factor as Jl but is of different magnitude. 
Usually Ri 0.4 to O. S although lower values are possible with 
5mal 1 baf f 1 e spaci ng . 

Sb = Correction factor for bypass flow CC and F streams!). It is 
iifferent in magnitude from Jb and ranges from O. S to O. S. 

=5 Correction factor for inlet and outlet sections having a 
raffle spacing different from central section. 

Before the correction factors can be determined it is 


lecessary to calculate various leakage and flow areas. 


These 
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auxiliary ealculatrioris ar^ bas^d upon t/bo ^oomo’Lry of* baf^f"lo 
gi ven i n appendi x C A3 - 

2. 2. 1 Shel 1 si de heat, t r ansf er eoeff'i ei ent 

To calculat^e t,he heat, trransf'er- ooef't^icient/ t,he st^ep 

procedure is as given below 

Ci3 Maximum mass flow velocity* G« 

* Amh 


* and 


wise 


C2. 23 


where Wi is total flow rate* and Amh is calculated and defined in 
Appendi x A* secti on C J 3 * equationC A, 1 73 . 

Ot Gim 

C23 Shell side cross flow Reynolds number* Re« = — ^ C2.33 

where Dt is tube diameter and is E)ynamic viscosity of shell 

fluid. 

C33 Ideal heat transfer coefficient for cross flow over bundle is 
given by 

hi = J - cp. G«. Pr« ^ C2. A3 

^here Pr«. is Prandtl number based on property values at average 
fluid temperature* cp is specific heat of shell fluid* is 

/iscosity correction factor* 

4,9 = ** C2 . S3 

/jw is dynamic viscosity of shell fluid at tube wall 
emperature. For a non viscous fluid such as water on shell side 
can be taken unity* 

is heat transfer factor defined as * J = St. Pr* ^ 
here St is Stanton numfaer and n is index of Prandtl number in 
the correlation of the type 

Nu = C. Re'". Pr'' 
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Value of hea-t -transfer factor j is reported in literature. For 

£ li. fit 3 

preaentr work it, has been t,aken from ehar-Ls prepared by 

Engineering Sciences Dat.a Unit^ Inlernat,ional Ltd. C ESDU > t and 
fitted in the form of logarithmic polynomials> 

In j = fCln ReZ) given in appendix B. 

C4D Having calculated Fe CEq A. 7Z> > appendix A» Segmental baffle 
window correction factor Je can be calculated from polynomial of 
the form» Je = fCFeD^ given in appendix B. 

CSD Using values calculated in sections ChD^ CiZ>» CjZ> of 



from polynomial of form* Jl = tCr im'J * given in Appendix B. 

r s 

Sirrdlarly baffle leakage correction for pressure drop* Rl is 


determined from polynomial Rl = fCrlwD^^* given in Appendix B 
C ©3 Usi ng Equat i ons CA.i4r3* CA.173* and C A. 1 1 3 f r om appendi x 


A* Following ratios are computed. 

Abe. 


C2. e:> 


N« = ^ C2. 

Nlee 

where Hm is number of pairs of sealing strips* shown in 
Cfig 4A3. Bundle bypass correction factor for heat transfer Jb 
is then calculated from the polynomial of the form* 

Jb = fCrb3.,‘*' given in Appendix B. 

Similarly bundle bypass correction factor for pressure drop Rb is 
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caleulated from polynomial of t.ypo* Rh = fCrb:>j^^^ » given in 

Appendix B. 

C7!> If shell side Reynolds number is less “than lOO* Cflow is 
laminarO then adverse temperature correction Jr is to be 
estimated- This is a function of the number of tube rows crossed 
as calculated* 

Ne = Wtee + Ntcw CS. 103 

where Ntee and Ntcv are calculated in Appendix A* and 


J V =*€ 


1 



for Re« > lOO 

for Re« <= 30 


C3. 113 


and a linear interpolation for 30 < Re« < lOO. 

CS3 Heat transfer correction for unequal end spacing * J« is 
deter mi ned f r om * 

. CNt - I'J + 

J« = C3. 133 

CNb - ID + Li'*' + Lo'*' 

where n = O. 6 for turbulent flow. 

For laminar flow correction factor Jc is about Cl+Js3y^. 
where L*i = Libl/L.b'O and L.o = Lbo/Lbe, 

Lbi. is inlet baffle spacing* Lbo is outlet spacing and Lbc is 
central baffle spacing as shown in fig.SA. 

The fig. BA. shows maximum unsupported tube span in central region 
is A* and in the inlet/outlet region it is B± and Bst. 

C93 Having determined all correction factors* shell side heat 
transfer coefficient is calculated by using equation 3,1. 


h« = hi- Je. J I- Jb- Jr- Js 


ca. 133 
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2.2.2 Stieil side pressure drop 

The shell side pressure drop from this met^hod is sum of 

cross flow pressure drop> APc, -the baffle window pressure drop 

APw ^ and fhe pressure drop in inlef and ouflef zones AP*. 

Ci3 Bundle pressure drop for ideal cross flow 

APbi. = 2. f . Mtcc. C2. 141) 

gc. 


where f is the cross flow friction f actor » gc is gravi tational 
constant Cfor SI unit gc = II). 

Cl2D Cross flow pressure dr op » 

APc = APbi.CNb - Rl C2-1B3 

where Nb is number of baffles ^ eq C A. 1 3D ^ Appendix A» and factors 
Rb and Rt are determined from polynomials listed in Appendix B 
C3D Window pressure drop » APw For turbulent flow CR*«i > lOOD 

APw = MbC2 O. S^WtcwDGw .Rt C2. iSD 

2. gc. 


where Gw is window mass velocity given as 

S- = 

jT I^Arrtb . Aw 

For 1 ami nar f 1 ow C Re « < 1 OOD » 

Gw . 


C2. 17D 


APw = Nb 


2S. 


pm 


r Ntcw ^ Lbc 

L» t p — Dt E)w J 


+ 2C10 i.Rl 


2ps 


C4D End zone pressure drop* 


AP* = APbi ( 1 + )Rb.Rs 


where 


r Lbc'i^-" ^ 

f Lbcl 

[ LboJ 

Lb V j 


Z-n 


CS. 183 


ca. 193 


ca. ao3 



where n = 0. 2 f‘or turbulent, flow and 1.0 for laminar flow. 
CB:> Stiell side pressure drop 

APs = APc + APv + APo C2. 21 Z> 

Pressure drop in inlet and outlet nozzles is not included here. 


2.3 


Kern method 


£123 


To begin with Bellas method » it is needed to have heat 
exchanger specif ications so that can be used for rating purpose. 
We can use Kern method to design heat exchanger > then this design 
is further checked and modified by Bell method » which is 
considered more accurate. 

2.3.1 Shell side pressure drop 

Stiell side pressure drop is proportional to number of times 
fluid crosses the bundle > for N number of baffles > 
number of crosses = H-*-l = Tube length /Baffle spacing > 
considering equal baffle spacing. 

The pressure drop on shell side can be written as 

f s . DfiC N +1 


AP® “ — 

where fa is shell side friction factor 
Gc shell side inass flow velocity 


C2. 


Gs 


Ac 


and» 


As 


Ds. CPitch-DO . Lbe 


Pitch 

Ds shell inner diameter 


ge gravitational constant, for S I tinit gc = 1 
ps shell side fluid density 



^0 


is correction f“or change in viscosity of fluid in 
houndary layer » For nonviscous fluids and turbulent flow 
= 1 » and E>* is hydr aul i c equi val ent di amet er for shel 1 si de 
flow. 

For this method the equivalent diameter D» is calculated 
considering the shell side flow parallel to the tubes. Tliough 
this method does not distinguish between relative percentage of 
right angle flow to axial flow. 


For QOt4rS deg layout angles 




4.CPt^- tiDtV43 
nDt 


For 3O>0O deg layout angles 


13« 


4,.CPT^it eos.30^ - nPt^y&J 
fiDt 


2S- 3. 2 Tube side pressure drop 

Tube side pressure drop is calculated by the following equation 


f t- Gt .L. 

2. ge- pi - 

where ft is tube side friction factor 


C2. 233 


Gt is tube side mass flow velocity 

2 


Gt = 


N t . n. Dt 


4. N tp 

Nt is Total number of tubes in exchanger 


C2. 243 


Ntp<s.«« is number of tube passes in the exchanger 
L is length of tubes to be used in exchanger 

<jfrt is tube side viscosity correction = 1 for nonviscous 

f lui ds- 


Shell and tube side friction factors are given by Kern as 
the function of shell and tube side Reynolds numbers Re«i and Ret 
respectively* which were fitted into polynomial for use in 



computer program. 
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Ree = 


Ret = 


CJa. & 9 , 

Dt. Gt 
P*- 


C2. 2S3 


C2. 26D 


where p* and pt are viscosity of shell and tube side fluids 
r espeeti vel y. 

2.3.3 Shell side heat transfer coefficient 

Tbe shell side heat transfer coefficient is calculated by 
the correlation suggested by Kern 

k« 




Jh. 


C2. 273 


Da 

where Jh is heait, *trarisf*er fac-ter given as a function of shell 

side Reynolds number > which has been fitted into polynomials for 

use in program» The Reynolds number used for estimation of Jh is 

same as used to calculate pressure drop. 

Pr« is Prandtl number of shell side fluid* Pr* = 

K •Sk 


3. 4 Tube side heat transfer coefficient 

For tube side also correlation given by Kern for water on 
tube side is used 

ht = Jhi. C2. SSD 

where Jht is heat transfer factor for tube side flow given as a 
function of Reynolds number *kt is thermal conductivity of tube 
side fluid* Prt is Prandtl number of tube side fluid and is 
viscosity correction factor for change of fluid viscosity* is 
taken = 1 for water. 



4 Overran Haa*t Transf’ar Coaf'f'iciantr 


Design Overall Heat Transfer Coef‘f*ieient lid is given as 


1 1 ^ Di 

4. Rfo + 

lid bi« 




1 .Dt 


C2. 293 


ht Di. ■ ■ ■ ■ 

where ha is hea-L -transfer- coefficient for outside of tubes,Ci.e. 
on shell side3» Rfo and Rfi. are fouling resistances on shell side 
and tube sides respectively* Ax is thickness of tube wall* Km is 
thermal conductivity of tube metal* Dt and EX are Outer and 
inner diameters of tube respectively. 


Clean Overall Heat Transfer Coefficient* Uc is given by following 
equati on 


1 ^ _1_ ^ 1 .Dt 

Ue ha ht EX 

neglecting the fouling factors on both sides* 

Total Fouling Resistance also called Dirt Factor Rd* 


C2. 303 


Ra 


Rfo Rfi. 


Dt 

EX 


so using equations Ca ZS to Z.'Sl 3 * 


C2. 313 


Rd = 


Ue - Ud 
Ue. Ud 


C2- 323 


2« 5 ENSfSlgn Considerations 

Tl-ie ESesign Equation is Q = Ud. A. Atm C2. 333 

where Q is amount of heat to be transferred* which is determined 
from the requirement of temperature change of fluids. 

Ud is E3esign overall heat transfer coefficient. 

A is total Heat transfer surface area* which is total outer 
surface area of all tubes* which is to be determined. 

Atm is mean temperature difference between shell and tube side 
fluids. For shell and tube heat exchangers with multiple tube 



patssfifs is deicer Kfti riSfd by f‘irs*t es'Llmait/irtg lo^gatrit/hmie 't^&mpf&ra.'ture? 
difference for a single pass counter flow exchanger and then 
using a correction factor Ft for the conf iguration under 
consideration. The value of Ft is always less than 1* given in 
the form of charts as a function of fluids temperature 
differences at ends and of the flow arr angements. 

The reason why logarithmic mean temperature difference is 
used and not the arithmetic mean is because arithmetic 
temperature difference is independent of the flow pattern 
both fluids flowing parallel or counter wise or across^ and 
remains same for all patterns. 

As shown by the design equation the compact design can be 
achieved by maximizing mean temperature difference and overall 
heat transfer coefficient. The counter flow arrangement gives the 
maximum mean temperature difference but this can not be used in 
all cases because of the following reasons- 

CID The tube lengths that can be used are limited by space 
available and handling difficulties. 

C2> similarly the maximum shell 'size is also limited because 
of cost of manuf acturing> transportation and reliability. 

C3I> Minimum tube side flow velocity is also limited to 
reduce fouling on tube side. 

so cost of ufianuf acturing and operating counter flow exchangers is 
the prohibiting factor, so are the shell and tube units with 
partially cross* parallel and counter flows in use. 

The usual procedure to get mean temperature difference is by 
the use of analytical procedures. Expressions for Ft are 
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availabl«& f’or 3 and 4 passes of tube in single shell but 

are limited in accuracy because of assumptions used for 
deriving^ like - 

Temperature at a cross section of shell and tube? exchanger 
is constant^ fluid property variation ignored^and numb&r of tubes 
in all passes is assumed samet are a few to list. 

A procedure has been written considering variation of number 
of tubes in different passes of a two pass shell and tube 
exchanger > in which fluid property variation is also possible. 
The method uses a piece wise integration of fundamental equation 
along the length of exchanger considering overall heat^ transfer 
coefficient constant along the lengthy thus yielding bhe overall 
heat transfer coefficient needed to achieve the given heat load C 
The amount of flow» and end temperatures of fluids I> with out 
actually calculating the mean temperature difference. The program 
has also the choice to calculate amount of area needed t,o achieve 
the specified heat transfer and overall heat transfer 
coef f i ci ent . 

Since the use of logarithmic mean temperature and correction 
for flow geometry^ Ft are used in design extensively and TEMA also 
provides the charts for Ft, it is also used as a check in program 
based upon Kern met-hod. It, is recommended t-hat, value of Ft 

should not, be less than 0.76 and also used as an criterion to use 
more number of shell passes than one. The reason is that values 
of Ft below 0.75 are in a steep region so with a slight 
operational transient it may change considerably to affect 
exchanger performance. Value of Ft increases with increase in 
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number of' ahell passes j. but, cost, and pressure drop also 
increases. 


Wlien 

■Lot-al 

area 

for heat exchanger » 

A is obtained 

in 

equation 

C2. 331) , 

total 

number 

of tubes hence the di ameter 

of 

shel 1 is 

deter mi ned for 

given 

tube di ameter 

» pitch* length 

and 

1 ayout angl e . The 

program uses 

data gi ven by 

Kern and Krause^ 

13 1 


f^or est,i matting shell diamet^er tctr given t,ube passes. 

The program TSHEET can also be used t,o calculate t,he Tull 
tube count for given shell diameter » tube diameter » tube pitch 
and 1 ayout angl e . 

TUBE CQUMT 

The number of tubes which can be accommodated within a given 
shell inside diameter is called tube count. It depends upon 
several design factors. 

Type of Exchanger : 

Floating head type exchangers have Outer Tube Limit COTLZ> 
smaller than fixed tube sheet type exchanger s» so tube count will 
be smal 1 er for same par amet er s C . e . pi t ch » tube di amet er etc . 3 . 
C2Z> Shell SJide Flow : 

Shell side flow determines the nozzle diameter and whether 
impingement plate is needed or not. Greater the nozzle diameter 
smaller the tube county so as to increase tube count externally 
fitted impingement plate or distributor should be used. 

C3Z> Number of Tube Passes: 

As Tube passes are increased area occupied by pass 
partitions increases so tube count decreases. 



C4D Tie rods. Spacers and Sealing devices also reduces ■the tube 
count . 

CS!> Tube End A't'taehmen't : 

Welded ends need larger margins, which decreases tube count. 
C&'J In some eases tube sheets may have radiussed hubs at the 
periphery Tor welding to the shell or head barrel, so as tubes 
not to encroach radii , OTL is reduced. 

C65 Rotatable bundles with internally fitted impingement plates 
requires two escape areas on opposite sides of the bundle. 

So it is not possible to provide exact tube count. Hie usual 
procedure is to estimate full tube count for given diameter and 
then considering above mentioned factors and design practice the 
full tube count is reduced appropriately. 

The Interactive Graphics Program TSHEET is prepared to calculate 
full tube count for given - 
tube sheet diameter 
margin from outer periphery 
tube outer di ameter 
tube SF»aeing or pitch 
1 ayout angl e , and 

location Ccoordinates^ of any one starting tube. 

The last variable makes the program very flexible that 
it is not necessary to have tube at the center , It can be forced 
to have tube at the location of designers choice with other 
variables remaining same. This feature is particularly important 
while designing for minimum leakage streams. 
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Fig. 6. THE ACTUAL ARRANGEMENT OF TUBES 
FOR 6 TUBE PASS. TEMA-S EXCHANGER, SHELL 
INSIDE DIAMETER 740 mm, TUBE OD 19.05 mm, 
PITCH RATIO 1 .33. LAYOUT ANGLE 30 





Fig.y. TOTAL TUBE COUNT 



siriee -there is no direct, procedure of* ge-t-ting actual tube 
count » the drawing of tube sheet layout with full tube count 
assists a lot in estimating of where tubes can not be installed^ 


thus arriving at the actual tube count. 

The Fig. 7. shows a computer produced drawing of 
layout with full tube count and fig. 6. shows 
arrangement for a 740 dia shelly 6 pass> S&O tubes » 
dia* for 2S. 4 pitch with 30 deg layout angle. 


a tube sheet 
the actual 
IQ. OB tube 



CHAPTER 3 


FLOW INDUCED VIBRATION 


3. 1 Reasons of Inclusion 

As a result of improved pressure drop predictability, 
allowing the use of higher fluid velocities the number of 
exchangers experiencing vibration are increasing. Higher fluid 
velocities and improper design are major causative factors. 

Flow induced vibration causes tubes to fall by following 
mechani sms- 

[13 Fatigue due to repeated bending. 

[23 ’Collision damage’ because of repeated impact between 
adjacent tubes; which leads to the thinning of tube wall. 

[33 ’Baffle damage* because of repeated impact between tubes and 
baffles, baffle tube clearance increases and tube wears out 
rapidly if baffle metal is harder than that of tube. 

[43 Cutting at the tube-hole edges at the inner tube sheet face 
due to repeated impact between tube and tube sheet. This also 
results in failure of tube-tube sheet Joint. 

The majority of failures have been found in regions of long 

and high local velocity, such as inlet 


unsuppor triSfd t/Ub<9' spains 
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baffle apace Ct,ubea in first, baffle window!?. Analysis of failures 
suggested bhat, had unsupported tube lengths been restricted to 
about 30$^ of TEMA values, few failures would have occurred. The 
effect of reducing the unsupported tube length to of its 

initial value increases the tube natural frequency by 56S4. 

3. 2 Tube Natural Frequency 

To predict the resonant conditions leading to damage it is 
desirable to estimate the tube natural frequency, which depends 
upon the type of supports and span length, is given as 

fn = 0.04,34,4 Cr. CE.I.gey^CW^ff.L'*^]*'® C3. 1!> 

where fn is natural frequency, E is modulus of elasticity of 
tube material, I sectional moment of inertia, L length of each 
span, Cn frequency constant Cdimensionless!? , effective mass 

per unit length of tube, ge gravitational constant C ge = 1 in SI 
unit Weff is the sum of mass per unit length of tube itself, 
mass per unit length of tube side fluid and virtual mass per unit 
length for the shell side fluid displaced by the tube. The 
virtual mass displaced is always more than actually displaced 
mass by a multiplier called ‘Added mass coefficient* CMci > 1>. 

The effect of axial stress upon natural frequency is not 
considered by above equation. 

For U bends the natural frequency is estimated by using 

longest bend length in above equation and then 

fnu = O. SS3 fr. C3. SD 

which is the frequency for out of plane vibrations. 

V.lu«= Cn -nrf H. .r. t^ker. 1“ =P‘''= 



3S 

eorist-ant, Cr. is almost, independent, of mode and =33. Ma is about 
2 to 3 depending upon pitch and layout of tubes. 

3. 3 Excitalion Mechanisms 

Several different theories are used to evaluate the 
potential for the vibration problem. These can be classified as 
follows- 

Cl assif ication of Flow Induced Vibration Mechanisms 

1 



Fluid Dynamic Excitation Motion Dependent Excitation 



Vortex Turbulent Acoustic Excitations 

Shedding Buffeting CWhen gas in shells 


Vibrations induced because of parallel flow have amplitude 
much smaller so not considered damaging. The mechanisms and 
methods important for liquid Csingle phased service are described 
under i n some detai 1 . 

3.3.1 Vortex Shedding 

When a bluff body is exposed to a perfectly uniform, 
nonturbulent incident flow above some critical velocity, the wake 
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behind the body develops instability* which oscillates in a 
regular sinusoidal pattern with an increase in velocity* and 
discrete vortices shed alternatively from the sides of the 
body*ealled * vortex shedding*. This produces alternating lift and 
drag forces on the tube. The frequency of the alternating force 
in the drag direction is twice that in the lift direction* but of 
smaller amplitude. The vortex shedding frequency in lift 
direction Cfv«I) is as given* by equation 3.3. When fv« approaches 
natural frequency it results in large amplitude vibration* called 
* Resonance * . 

f V<ft — iSlc. Vc/EJo Cd. 

where Sle = Strouhal number* Ve = cross flow velocity based 
on minimum transverse area between tubes in a row* E>o = tube 
outside diameter. 

When tubes vibrate at* or near* vortex shedding frequency* 
vortices shed at the tube natural frequency which is not given by 
the above equation. This feature is called lock ""in effect* 
which can prevail at a frequency of about 0. 7-1 . 3 of the tube 
natural frequency. Therefore design criterion is to keep 

fn/fva ^ O. S. 


3. 3. a Tur bul ent Buff ebi ng 

Unsteady forces develop on a body exposed to highly 
turbulent flows, having a wide band of frequencies around a 
dominant central frequency which increases with the flow rate. 

No data or equations are available for turbulent buffeting 


frequencies for liquids. 
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3.3.3 Fluid elast,ie whirling 

It, is a tube motion induced excitation mechanism. A 
momentary displacement of" one tube in an array I’rom its normal 
position alters the flow field, this disturbs the force balance 
on neighboring tubes and exposes them to a fluctuating pressure 
and causes them to change their positions in a vibratory manner 
at their natural frequencies. Thus this flow induced vibration 
originates from fluid elastic coupling. Tube motion may begin by 
vortex shedding or turbulent buffeting and become a catalyst for 
fluid elastic vibration. 

When fluid velocity is less than a critical velocity 
excitation forces are balanced by fluid damping. Critical cross 
flow velocity can be calculated by the empirical equation 3. 4, 

Vent = ft. f ri. I!>e.. y ^Wof f ■ <So C3. 4!^ 

Verii is crit-ical cross riow voloei'ty* ft ins'tabi li*ty 
corist.anL* 6o logaritrhmic decrement, > Da t.ube OD and p shell f’luid 
densit,y. Sa = 0.1 IT or liquids^ and 0.03 f*or gases » W*ff Is 
defined in section 3.2. 

3a 4 Damage Numbers 

Thor ngr en pr- oposed two di mensi onl ess number s to be used 
to predict damage due to cutting action at the baffle or to 
collision between tubes . 

Do. Vci 




Fb. Sf- ge. Am. Bt 


C3. S3 
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Ncb = 


Fb AwC Do* 4-EX^!!>C t . E. ge 


C3. 03 


Hbb is fcja.f*fle cianfta^e* riumtj^Sfr B.rid Nob is collision cia.i]^go 


number. where Fb 


bube *to baff’le clearance factor » 


di mensi onless. Fb — 1.00 for 1,«^3 in* clearancet 1 . SB for 1/64 in* 
clearance. Sf = fatigue stress > Ct = minimum gap between adjacent 
tubes |. Am = cross sectional area of tube metal * Do and EX are 


outer and inner diameter of tube> Ve = cross flow velocity^ p« is 
shell fluid density* L = tube span length and go is gravitational 
constant C= 1 in SI unit3. 


The damage numbers are indicative of damage which vibration 
might cause* and Thorgren** ^ stated that this would be avoided 
if both damage numbers are less than unity. 

C B 1 

ErskirnOf and Waddingt-on modlf"l«rd -thef abov^ design enit-erion 

by studying case histories, that damage would be avoided if: 

O. 004,20 p^' 

Nbd and Won < C3.7I> 

O. 235 

/J« 


For vibration calculations cross flow areas are based on the 
free area between tubes in same row. The cross flow areas are 
therefore identical for those for heat transfer and pressure loss 
calculations for 30 and 0O pitch angles, but differs for 4B and 
60 pitch angles. 


3-5 Etesign procedure 

The step wise procedure used by the computer program VIBUNIT 
to predict the possibility of damage because of flow induced 
vibration is as described under. 



1. Ca.lcul 1,11© t,ube naLural f^r^qu^ney by using ©qtiabion 3.1 
tctr longest, unsuppor'ted span of’ t,ub©. . The longest, unsupport,ed 
spans for inlet, > out^lel and central regions are Bi» B*, A as 
shown in figure 5 A. 

3. Caleulat,e vor'tex shedding frequency by equation 3.3. 

3. If ratio of fv« to fn is less than O. B ^ the design is safe 
from the vortex shedding vibrations^ otherwise corrective 
measures include the guide lines given in next section to make 
the tubes stiff or to reduce shell flow velocity. 

4. Calculate critical cross flow velocity given by equation 3.4. 
If maximum cross flow velocity is smaller than critical cross- 
flow velocity the design is safe. 

5. Calculate Baffle* Collision and maximum damage numbers from 
equations 3. B to 3.7. For safe design Baffle and Collision 
damage numbers should be smaller than maximum damage number or 1 
whichever is smaller. 

3« 6 E)esign Recommendations for vibration prevention 
Flow induced vibrations can be prevented by 
Cl 3 Increasing tube natural frequency* and 
C i i 3 Reducing cr oss f 1 ow vel oci t y . 

Following are the recommendations- 

[13 A cross flow shellC X type 3 or a no tube in window design 
virtually guarantees freedom from vibration failure. 

[23 If damage is predicted near first row of tubes at inlet the 
problem can be solved by installing a support plate in the inlet 


baffle space. 
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133 Shell side flow velocity can be reduced by using double or 
triple segmental baffles or by increasing shell side passes or 
by the use of wider pitch and pitch angle. 

[43 U bend can be supported or all tubes at U bend can be clamped 
together to give more stiff design. 

[B3 Reduce baffle tube clearance and use baffle material softer 
than that of tubes ^ and increase baffle thickness. 

[63 Mew design concepts such as Rod baffles or MElST supports can 
be considered. 



CHAPTER A 


PROGRAM DESCRIPTION 


The Pa.sce.1 programming language has been used in all 
implement, at, ions at t,his project,. There are t,wo execubable f'iles» 
namely HED. EXE and TSHEET. EXE. This chapter describes the uses* 
available features and limitations of both files. 

The following four routines are available from the main menu 
driven program, HED. EXE: 

C1I> Rating : 

This routine can be used to calculate overall heat transfer 
coefficient required for an exchanger for given end temperatures 
and surface area of heat exchange. The procedure uses piece 
wise integration over the entire surface of heat exchanger to 
calculate the overall heat transfer coefficient. This makes it 
possible without calculating the logarithmic mean temperature. 

This program can also be used to calculate required surface 
area for an exchanger for given end temperatures and overall heat 
transfer coefficient. 

Tlie major assumption in above formulation is that the 
temperature changes only along the length of heat exchanger and 
not across it. In context of cross flow in a shell and tube heat 
exchanger usually with more than 10 baffles the assumption is 
appropriate. 
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This routine has two specialties^ 1. It can handle variation 
in specific heat of fluids with temperature. 2. Usually numher of 
tubes in different passes of tubes differ by 5 to lO ?S and are 
not ec^ual . This routine can handle different number of tubes in 
different passes. 

The routine is limited only to the 2 tube pass or a U tube 
shel 1 and tube exchanger . 

C25 Kernes Method : 

This routine is principally used in our calculation to 
arrive at a first design. The inputs are the fluid properties » 
tube parameters » and number of shell and tube passes. The 
complete formulation is given in chapter 2. This routine has its 
database contained in files named KERN.?OOf where stands for 
positive number. File KERN. lOO is containing data to convert BWG 
thickness of tube to thickness in inches. Files KERN. 200 to 
KERN. 700 contain data to get the shell diameter and tube count 
values for different tube diameters and tube patterns. Ctoviously 
the data is limited and can be added if available for more tube 
diameters and patterns. 

C3> Bell Method s 

Tlie routine is available independently as well as through 
the Kernes Method. 

For the purpose of rating of an existing heat exchanger it 
is necessary to call it Independently from the main menu. It 
permits data input from a file as well as from screen. To begin 
data can be read from screen and when input from screen completes 
it makes a file of input data, BELLI N which can be used next tine- 
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if* more runs a.re required. 

For designing & heaft exchainger it^ should be ea.lled f‘rom 
Kernes Method after executing Kernes method. This generates a 
preliminary design and this will then be passed to Bell method 
for further modif ications. 

The formulation of the method is given in chapter 2. 

The method contains the data regarding the geometrical 
clearances of four type of heat exchangers for TEMA class R 
exchangers » so it can only be used for these four types of 
exchangers - 

1 Fixed tube sheet type heat exchanger 

2 U tube heat exchanger 

3 Split backing ring type floating head heat exchanger 

4 Pull through type heat exchangers. 

The various types of heat exchangers are introduced in 
chapter 1 with their principal features which will help in their 
selection while running the program. 

C4> Flow Induced Vibrations : 

Tills unit is also available independently as well as through 
t he Bel 1 met hod . 

When it is needed only to perform vibration check for an 
exchanger than this unit can be called independently* otherwise 
the call from Bell method will generate the vibration check list 
for the design under consideration. 

Tlie complete formulation is given in chapter 3. This unit 
checks for the possible damage by vortex shedding* fluid elastic 
whirling* and for the collision and baffle damage numbers. If 
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shell side fluid is a gas -than -this will also cheek for* *the 
acoustic vibrations additionally. 

The file TSHEET.EXE can be used to calculate the full tube 
count. The actual tube count differs from the full tube count 
because of the reasons given in chapter section S. 4. 

The inputs are tube and tube sheet diameter ^ tube pitch and 
layout angle. The margin from outer edge can be changed. To 
enhance the creativity the coordinates for the tube to begin can 
also be altered. This feature might help in obtaining final 
tube sheet layout design for rrdnimum leakage through the pass 
partition lanes in a heat exchanger with more than one tube 
passes . 

This program is graphics implemented so the designer can get 
drawing of the actual tube layout ► by deducting the number of 
tubes which has to be removed to provide pass partitions. 

Using this software » drawings can be stored * viewed and 
compared together . 



CHAPTER S 


DESIGN PROBLEM AND RESULTS 


"THis cHapt.er ciescribi&s t/Ke ^design crorisidera^ci 1,0 

analyze. The known parameters of the design are as follows : 

Hot stream consists of water to be cooled from BS *C to 4B 
C. The cold stream is cooling water available at 3B *C » which 
could be heated up to 40 C. The mass flow rate of hot stream is 
25 kg/see. Length of tubes in exchanger is limited to 4 meter 
because of the limited space available for handling. 

The thermodynamic properties of fluids used are as follows: 
Specific heat for water on shell and tube side = 41 S6 J,^kg. K 
Tliermal conductivity of water for shell and tube side = 0.62 W^^mK 
Dynamic viscosity of water on shell side = S.l E - 4 N.sec^? 
Dynamic viscosity of water on tube side = 9.1 E - 4 N. sec/^m*. 
Specific gravities of water on both side = 1. 

With the above information the first run of Kern‘s method 
has taken considering hot stream in shell side. The reason to 
assign shell side to the hot stream and the tube side to the 
cooling water is mentioned in chapter 1. The layout of tubes is 
taken triangular considering the hot water clean. The tube 
diameter 0.75 in.Ci9.05 mwO and pitch 0.9375 in. are selected 
from the available database of the program» because of the 
standardization reasons. The tube count and shell diameter is 
then calculated and taken from the nearest available standard 
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diameter. Knowing shell diameter CB40 rmO the baffle spacing is 
selected C3BO mm3, as suggested by baffle spacing can be 

from one fifth of shell diameter to shell diameter. This 
completes the calculations of Kern method. The data calculated 
is then transferred to Bell method and thereafter to for the flow 
induced vibration routine. 

The effect of varying the various shell side parameters on 
the shell side pressure drop and heat transfer coefficient is 
shown in plots Cfig.6 to 313 for the heat exchanger under 
consideration. These plots can be used for the selection of 
baffle spacing , baffle cut and number of pairs of sealing strips 
for optimal heat transfer coefficient and pressure drop. 

EeXection of baffle spacing s 

Fig. 6. shows that the shell side pressure drop is 
decreasing exponentially for increase in baffle spacing from one 
fifth of shell diameter to that equal to shell diameter. Ttie 
plot shows that decrease in shell side pressure drop is not 
significant as the baffle spacing is increased from 250 mm. 

The fig.Q. shows the corresponding decrease in shell side 
heat transfer coefficient by increasing baffle spacing. The 

diagram shows a sudden rise in heat transfer coefficient for 
baffle spacing about 365 mm. The reason is that after a baffle 
spacing of 365 mm, according to TEMA recommendations the baffle 
hole to tube diameter clearance is reduced to from l.-'''32**. 

This results in reduction of leakage stream A Cfig. 4.3. In our 
computation the value of factor Jt. increases from about 0.72 to 



84 


O. SI » hence the sudden rise in shell side heat transfer coefficient 

Both of the above discussed fig. S and fig. 9. are plotted 
considering there is no sealing strips * i.e. Hm = O. 

Selection of Baffle Cut s 

Fig. lO. shows the decrease in shell side pressure drop with 
increase in baffle cut. The drop is proportional , 

Fig. 11. shows the variation of shell side heat transfer 
coefficient for increase in baffle cut. T^ie graph shows that 
drop in shell side heat transfer coefficient is more steep as 
baffle cut is increased after 20 of shell diameter. The reason 
is that for small baffle cut Cwith small baffle spaei ngZ> the 
amount of leakage streams is higher and as baffle cut and spacing 
are increased the leakage streams also reduces because of 
reduction in flow resistance. As baffle cut and spacing are 
further increased the amount of parallel flow increases which 
also reduces the shell side heat transfer coefficient. Thus the 
baffle cut of about 20 % of shell diameter seems appropriate 

particularly for small baffle spacings. For larger baffle 
spacings 205^ baffle cut can be used. 

Selection of Number of pairs of sealing strips s 

Fig. 12. shows the drop in shell side pressure drop with 
increase in baffle spacing considering two pair of sealing 
strips. Similarly fig. 13 and fig. 14 shows variation of shell 
side heat transfer coefficient for increase in baffle spacing and 
baffle cut respectively considering two pairs of sealing strips. 
The trend is almost same as it was without sealing strips in 
fig.S to fig. 11. The fig. IB shows the decrease in pressure drop 



wi'th incr«&a.so in ba.f*f^le cu*t> As «&xp©et,erd t.h^ d€Fcrea.s& is s'Leep 
for smaller baffle spaeings Cfig.iB:^, than for the larger baffle 
spaeingsCf ig. 163 . The curve shows a comparati vely steep fall in 
pressure drop for baffle cut more than AO ?£. For s*uch larger 
cuts the shell side flow essentially becomes parallel that is why 
the drop in pressure drop. The shell side heat transfer 
coefficient also suffers a steep fall for large baffle cuts 
because of same reason Cfig.lA3. 

Fig. 17. shows the variation of shell side heat transfer 
coefficient for increase in sealing strips. As sealing strips 
are provided to reduce the amount of bundle bypass streams. The 
curve shows that the rise in heat transfer coefficient is 
insignif leant for more than 5 pairs of sealing strips » and not as 
steep for more than 2 pairs of sealing strips. 

Fig. 19. shows the rise in shell side pressure drop for 
increase in sealing strips > the trend is similar to that of the 
heat transfer coefficient Cflg.173. 

Designing for compactness : 

Fig. 16 to 21 shows the variation of full tube count for 
variation in shell diameter for different tube diameter* pitch 
ratio and layout angle respectively. These curve are produced by 
data generated from Interactive Graphics Program TSHEET. These 
curves depicts that the compactness C larger tube counts* is 
affected maximum by the change in tube diameter and change in 
layout angle affects minimum. So as far as the selection of 
physical parameters is concerned the change in layout angle must 
be tried first* change in pitch ratio next and change in tube 



diameter the last. 


The output specifications of the designed heat exchanger is 
given in Appendix as is obtained from the computer. It 

contains output from the three routines » The Kern*. Bell» and 
finally the vibration check report from Vibunit. 



CHAPTER & 


CONCLUSIONS AND SCOPE 


An i n'tG'r Act.i Vifr > m^nu driven code f*or t/hermaLl hiydra.ulic 
design end reling shell end ■tube hee*t exchengers hes been 

developed. The epplicetion "to t,he typicel design problem 

produced results were logically discussed in previous chapter. 
Considering the importance of flow induced vibration check ing» 
right et the time of thermal hydraulic desi gn» the routine has 
been added » which is also accessible independently from rest of 
the code, 

Thie pr ogr am is sui t abl e onl y for si ngl e phase exchanger s , 
with plain tube only. The liquid metal flow also not considered. 

Inclusion of correlations and data for low fin tubes and two 
phase flow as well as for liquid metal heat transfer can further 
enhance the utility of this software. 

Another improvement could be addition of optimizing routines 
to let computer make decisions about optimal parameters like flow 
velocity^ cooling water temperatures ► and area for heat transfer 
etc. Yet another possibility is of interfacing mechanical design 
and cost analysis codes and globally optimizing the design. 

Ref ’^"^^and^'^^ also suggest additional features to supplement 


a design code. 
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APPENDIX A 


AUXILIARY CALCULATIONS 


Ca.I> Segmen’tal BafTle Window Cal eul a't ions 

The central angles of* the baffle cut with the inside of 


shelly ^dc. and the angle intersecting the tube diameter 

Ektl = Dotl - Dt are C see fig. 22 Z? CA. 15 

= 2cos Cl— 2Bc/''i005 CA. i:35 

= a=«-‘ { 1 )] c A. 33 


where Be is baffle cut expressed as percentage of shell inner 
di ameter . D« i s shel 1 i nner di ameter > E>otl i s outer tube 1 i mi t as 
shown geometrically in fig,S2. 

If baffles are of window type»C generally used in floating head 
type exchangers!) then area between Dc and Efetl is blocked » so the 
angle is referred to 

doll = )] CA. 45 

Cb3 Baffle Window Flow Areas 


The gross window flow area» i.e. without tubes is 

. n - 2 : sin ^ i 

^''3 = — D' J CA.aJ 

If window baffles are used t/hen Sda in above equation is replaced 
by &otl. Assuming a uniform tube field the fraction of tubes in 

one window »Fw and the fraction in cross flow between baffle tips 

_ dell sin dell CA. 65 

Fc are Fw = ^ 

Fc = 1 -2Fw C A. 75 

If tube field is not uniform because of pass lanes or tubes 
removed in the nozzle entry zone then this correction has to be 
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sseparatrely calculated. 


The baffle window area occupied by tubes* Avt ii 


Awt = Nt Pv 


nDt 


CA. 23 


and the net flow area through one window* Av is 

Av = Awg - Avt CA. 93 

where Wt is total number of tubes in exchanger* 

Dt is tube outside diameter. 

Cc3 Equivalent Hydraulic Diameter for the Window* Dv 

This value is required only for pressure drop calculation in 
laminar region CR* < 1003. 

rt. ^ V * ^ 

n Dt Fw Nt 4. n D^rf:aa«x2fz3 CA. 103 

neglecting the baffle edge thickness. 

Cd3 Number of Effective Tube Rows in Cross Flow* Ntcc The number 
of tube rows is a function of tube layout and pitch parameters. 
Lpp is defined in fig. 23. For one cross flow section between 
baffle tips 

= S? ( ^ ■ rfe ) 

C&J> Number- of" ve Tube Rovms in Window, Mlcv 

It, is basod on int.orpret,at,ion oiT the window fTow pattern, 
data Tr-om E!>el aware method can be approximated as 0.4 times the 
tube field in the window. This distance is crossed twice in each 
wi ndow V hence 


Ntc\ 


0 . 8 
Lpp 


Db 


Be _ Db - Detf 
lOO 3Db 


CA. 12'.> 


C f Z> Number of Baf f 1 es , Nb 


For calculation purpose number of baffles is determined from 
tube length and central baffle spacing even if larger end 
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spacings are used. The ■Lube lengt-h Lii is defined in fig. SB. 

^ CA.135 

Cg'J Bundle Shell Bypass 

The bypass area within one baffle Abet* is 

Abet = LbeC CD«~E>oiD+ Lf>l3 CA. i4D 

where Lpl is the tube lane partition bypass width as follows: 


O for all standard calculations. 

Lpt = < 1/2 of real dimension of tube lane partition* Lp. 

^ For estimation purpose Lp = Dt. 

C h3 Shel 1 —Baf fie Leak age Ar ea » A«b 

The diametral clearance between the shell diameter D« and 

the baffle diameter Db is designated Lsb. The shell baffle 


leakage area within the circle segment of the baffle is 

L s dfT 


A«b = fiDa 




CA. iB'J 


where L«b = Da - Db 
Cil) Tube Baffle Leakage Area* Atb 

The diametral clearance Ltb is determined from specified 
manuf acturing tolerances* which are given in TEMA. Ttie total 
leakage area for a baffle is 

Atb = ^ ECDt + Ltb'J®- Di®3NtCl-Fv:) CA. 163 

A 

CjZ) Minimum Bundle Cross Flow Area, Awb 

Tills is the minimum cross flow area for one baffle at the 
shell center line determined as follows: 

Amb = Lbcjhbb + - Doj CA .173 

where Lbc is the central baffle spacing* CLtb^aff is Ltp for 30 
and 90 degree layouts * 0. 707Ltp for 45 deg. layout and O. SSSLip 
for 60 degree 1 ayout * C shown in fig. 23. 3. 



APPENDIX B 


Thie appendix liats the polynomial a used by r-outine BELLUNIT 
► the terms used are defined in chapter 2. The polynomials are 
fitted by best fit method for specified range only. Function 
F>owCx»nZ> is used to represent x". 

Friction Factor f and Heat transfer factor J for shell side 
Calculations 

For layout angle ■ 30 deg 

Res >= iO and Res <= 300 logC j3 : =3. 57a7e-l +3. 6SB2e-i4elogCReO 
Res > 300 and Res <= leS logC jU : =-l . 34S3+6. 3731e-l»logCRe^ 
for pitch ratio = 1.25 

Res >=10 and Res <=300 logCf:j = 1.3311 - 1.5163 « logCResO + 

8. 0455e-2 « aqr Cl ogC ResI) :> 

Res >300 and Res <= ie4 iogCf> : =2. 377 - 1.0013 * logCRes> + 

5. 0484e-2 « sqrC logCRes>> 

Res > le4 and Res <= le5 logCfZ) = —2.2246 + 2. S706e— 1 * logCRes> 
-2.SllSe- 2»faqrC logCRes>> 

Res > le5 and Res <= le6 logCfD = 3.7237 - 8. S957e-1 * logCResi + 
2. 3240e-2 » sqrClogCRes>> 

For pitch ratio = 1.33 

Res >= lO and Res <= 300 logCf> = 4.6744 — 1.5446 w logCRes!) + 

9. 1024e -2 « sqr Cl ogC Res> > 

Res > 3CX5 and Res <= le4 logCf> = 2.4338 - 8. 9004e-l » logCRes> + 
4.38S6e-2 * sqrClogCRes>> 

Res > le4 and Res <= le5 logCf> = 1.2393 - 3.3821e-l a logCRes> + 


4. 579Be-3 » sqrClogCRes>> 

Res > le5 and Res <= le6 logCf:) = 1.3153el- 2.444 « logCRes3 


4. 


7S 


9. 2416e-2 » sqr-ClcigCKe05 
For pitch ratio = l.S 

Res >= 10 and Res <= 300 logCf^ = 3. B34 - 1.3723 « logCResO 

+S. 7BSSe- 2 « sqrClogCReOi 

Res > 300 and Res <= le4 logCfD = 2.3054 - 1.0729 « logCReslJ + 

5 . S273e— 2 «sqr C 1 ogC ResD Z) 

Res > le4 and Res <= le5 logCf:^ = -9. 0367e-l - 1 . 5826e-2 » 

logCResD - 1 . 4204e-2 » sqr ClogCRes35 

Res > le5 and Res <= leS logCf5 = -2. 367e-l - 3. S424e-1 * 

logCResD + 1 . 0979e-2 *sqr ClogCResD5 
For layout angle ■ 45 

Res >= lO and Res <= 300 lnCJZ> = 2. S37Se-l + 3. 5609e-l » logCRes5 
Res > 300 and Res <= le6 lnCJ3 = -1.4142 + 0. 4911e-l * logCRes^ 
for pitch ratio = 1.25 

Res >= lO and Res <= 300 logCfi = 3.9253 - 9. 43e0e-l k logCRes^ 

Res > 300 and Res <= leS logCf5 = -S.754e-1 + 1 . 2664e-l » 

logCResD — 5. 3203e— 2 i<sqr ClogtResZ>Z> + 2. 3701 e— 3 w powC 1 ogC ResZ? » 3!> 
for pitch ratio = 1.33 

Res >= 10 and Res <= 300 logCfD = 4.0674 - 1.2915 * logCResZ) + 
5. 49Se-2 « sqrClogCResDS 

Res > 300 and Res <= leS logCf3 = -4.4619e-l - 6. 03e-2 » logCReslJ 

- 3. 2313e-2 « sqr ClogCResD^ + 1 . 67S5e-3 * powC 1 ogC ResZ) » S!) 
for pitch ratio = 1.5 

Res )► = 10 and Res < = 300 1 ogC f Z> = 3 . 4945 — 1 . 3659 w 1 ogC ResZ> + 

8. 2815e-2 « sqrClogCResDZ) 

Res > 300 and Res <= le6 logCfS = -1.7959 + 2. 4112e-l « logCResZ) 

- 5. 37S9e-2 sqrClogCRes^D + 2.1602e-3 » powC 1 ogC Res5 » 35 
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For layout, angle ■ QO 

Kes >=10 and Re« <= 300 lnCJ> = -3. 064e-l + 4. 2S53e-l * logCEe=> 

Ke« > 300 and Res <= le6 InCJ'J = -1 . 636B + 6.SQ61e-l » logCRes> 

tar pitch ratio = 1 . 3S 

Res >= lO and Res <= 300 logCf> = 4.0S16 - Q.QlOSe-l » logCResD 
Res > 300 and Res <= 3000 logCfJ = S. 64i3el - 2. 3691 el * iogCResS 
+ 3.234S « sqrClogCRes>:> -1 . 454e-i » powCl ogC Res> » 3’J 

Res > 2000 and Res <= 2e4 logCfO = 1.01 S2el - 4.682 » logCRes> + 
6. 3149e-l « sqrClogCRes>> - 3 . ee22e-2«powClogCRes> » 3> 

Res > 2e4 and Res <= leS logCfO = -1 . 2899el + 3.3987 « logCRes> - 

3.3493 e-1 w aqrCiogCRes>> + 1 . 0271e-2 « povrf: 1 ogC Res> » 3> 
f'or pitch ratio = 1.33 

Res >= lO and Res <= 300 logCf> = 3.5777 - 9. S343e-1 » logCRes> 

Res > 300 and Re= <= 2000 logCf> = 3. 904eel - 1 . 73Bel » logCRes> + 
2.4138 « sqr ClogCRes» - 1.113e-l « powC 1 ogC Res> » 3> 

Res > 2000 and Res <= 2e4 logCf> = -1 . 5003e-l - 1.3218 » logCRes> 

+ 2. 2156e-l « sqrClogCRes>> - 1.1061e-2 « powC 1 ogC Res> . 3> 

Res > 2e4 and Res <= leS logCfl = 5.6469 - 1.2588 « logCRes> + 

4. 6627e-2 »* sqr ClogCRes>> 
for pitch ratio =1.5 

Res >= lO and Res <= 300 logCf> = 2.8302 - 9. 4005e-l » logCRes> 

Res > 300 and Res <= 2000 iogCf> = 4.0023el - 1 . S296el * logCRes> 

+ 2.5919 » aqrClogCRes» - 1.2107e-l » powC 1 ogC ResD , 35 
Res > 2000 and Res <= 2e4 iogCf> = -9.8516 + 1.6407 * logCRes> 
7.9541e-2 » sqrClogCRes» - 1.1427e-3 « powC 1 ogC Res> » 3:> 

Res > 2e4 and Res <= le6 then logCO = 4.1080-1.0597 « logCRes> 
3.9501e-2 « aqr C 1 ogC Res> > 
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f = expC 1 ogC fO D 
Jh = expClogCJ^'J 
j = J h/'Rea 

Factor-s Jl, Jb, Ri, Rbf, and Jc for- shell side heat, t.ransf'er 
coefficient, determi nation 
for rim >■ 0.0 and rim <■ O. S 
for rs = O 

Jl = 9. 917e-l - 1.117p«rlm + 9. a394,e-l «r-lTO»rlTO - 2. 9242e-i « 

powC r Iwi 1 3^ 

for rs = 0. 25 

Jl = l.OOei - 1.49134erlm + 1 . 48S1 «r ImXr Im - 6. 2222e-l » powCrltn.S^ 
for rs = 0. 5 

Jl = 9. 9736e-l - 1.6064«rtTO + 1 . 431 »r lm»tr Im - 5. 2773e-i « 

powC r im » 3Z) 

for rs = 0. 75 

Jl = 1.01S6 — 1.944lwrlm + 1 . 731 3*r li«wr Iw — 5. 7273e— 1 w powCrtn>»3I> 
f or r s = 1.0 

Jl = 9. 92Sle-l - 2. 1404itrlm + 2. 124S«r lm«r Im - S. 7955e-l « 

powCr lm,33 

for intermediate values of rs linear interpolation of upper and 
lower values is used. 

factor Rl 

for r I to >= 0.0 and rim <= O. S 
for r« = 0» 


rl 


1 . 0235 


4.3203 * rim + 2. 0627el » rim * rim 


6. 3964el * 
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powCrtm ,33 + l.lisaea a powCr'lm. 4:5 - l.OOSlea * powCrlY«,B3 + 

3 . 61 a4el » powC r Im » 63 

for rs = O. as , 

Rl = 8. S113e-i - a. a466 » r Im + 4.3365 * rtm m rtm - 4.8543 » 

powCrlm»33 + a. 1875 « powCrlm»43 

f’or r« = 0.50 » 

Rl = 7. 581 4e-l -1 . 9554 k rtm + a. 5939 a rim « rim - 1.3358 * 

powC r Im » 33 

for rs <= O. 75 , 

Rl = 1.0081-9.1831 « rim + 6.301601 » rim « rlm-a. 491503 » 
powC rtm, 33 + 5. 8633ea » powC rim, 43 - 7.961703 w powC rim, 53 + 

5. 770303 4* powC rim, 63 — 1 . 73S7o3 » powCrlm,73 

for r « = 1.0 , 

Rl = 6 . OI 60 -I - 1.9411 » rim + 3.8339 « rim « rim - 1.4813 « 

powC r im , 33 

polynomials for ra = 0.35 and ra = 0.5 aro not acourato NEAR rim 
= O, but aro aceurato for rim = 0.1 to O. S so for r-lm = 0.0, Rl = 
1.0 is usod soparately. for int or modi ate values of rs linear 
interpolation of upper and lower values is usod. 
factor Jt 

For Slioll side Reynolds number > 100 and rb >= 0.0 and rb <= 0.7 
for Naa-f >= 0.5 Jb = 1.0 

for Naat = 1.0^3. O Jb = 1 . 0031 -1 . 4386e-l » rb 
for Naa+ =0.3 Jb = 1 . 0014-3. 9643e-l * rb 

for Naa+ = 0.1 Jb = 1.0107 - 5. 39SSe-l » rb + 1.3690-1 » rb « rb 
for Naa+ = 0.05 Jb = 1.005 - 6. 6439e-l » rb + 1 . 7857e-l « rb « rb 
for Naa+ =0.0 Jb = 1.0036-1.197 * rb + 5.39760-1 » rb « rb 



for i ritrormodi a.tr<& valuos of* rs linear int.erpola'Liori of' upper and 
lower values is used, 
rb = O. O Jb = 1 . O 

Kecomiifiended lower limib of* Jb is given fcjy f'ol lowing polynoiftial 
Jb limit = ©. 3Be-i + 1 . 0734e-i m rb + 4.©©le-i ^ rh ^ r b-i . 3©8Qe--i 

w powC r b » 3D 


f“aetror rbf 

Applioabili ty linfdt rb >= 0.0 and rb <= 0.7 
f'or N»»+ >= 0,6 Rbf = 1.0 

for Nfiftf = 1. 0/3.0 Rbf = 1 . OOSS - 4. 9464e-i m rb + 1 . 36e-l ^ r-h ^ 

rb 

for N«ish. = 0.3 Rbf = 9. Q7S6e-l-9- aaoae-l ** rb + 3. i64Se-l m rb » 

rb 

for Nc©+ =0.1 Rbf = 1.000 - 1.4013 m rb + 6. 307ie-i it rb m rb 
for = 0.05 Rbf = 9. 93e~l ~ 1.3014 m rb + 1.0714 it rb it rb 

for Nw =0.0 Rbf = 9. S6e-i - 3.176 it rb + 3. a6 it rb it rb 
for rb= O rb f = i 


factor Jc 

je = 5. 3539e-l + 3.6331e-l it fe + 4.54630-1 it fc it fe - 4. S66Se"l 


it powCfG»3D 



APPENDIX C 


SPECIFICATIONS FROM KERN METHOD JM* 


‘.ONSTRUCTION DETAILS: 


Size 


7S C HmO sc 4 . O m 


Number of shell passes 
Number of tube passes 


Number of Tubes 


Tube Outer Diameter 


IQ. OB Cmml) 


Tube Inside Diameter 


le.seoe cmmu cis bws3 


Layout Angl e 


30 deg. 


Pitch 


23. S12S CmnO 


Tube Length 


Baffle pitch 


2SO mm 


Number of Baf f 1 es 


Tliermal conductivity of tube metal 


112 WXm. K 


THERMO HYDRAULIC DETAILS 


SHELL SIDE 


TUBE SIDE 


uid Circulated 


C Kgysec'J 


eeific Gravity 


scosi ty 


C N . sec/^sq_HD & . 1 E-4 


Q.lE-4 


eeific heat 


CJ/Kg. K5 


4.186 


4186 


ermal Conductivity KZ) 


0.62 


O. 62 


Kiperature in 


'emper- atur e out, 


CO 


45 


40 


CLlowable pr-oss.dr-op CN/sq_KO 1 . OE5 

:aleulat,ed Press, dr- op CN,-'aq_HD S. 225E4 
‘ilm coefficient CW>^sq_m. K5 6215 

:ieari overall heat transfer coefficient CW,-'sq_m. K5 
iesign overall heat transfer coefficient CW,-'sq_m. K> 
:al cul at ed C avai 1 abl e3 Di r t Factor C sq_m. K„^W5 

ube side velocity Cm,^sec. 


1. OE5 

2. i23E4 
5403 

= 2664. 7 
= ilOO 
= 5. 335E-4 
= 1 . 3558 


C BELL - TABOEEK METHOD 5 
Additional Specifications 
U TUBE TYPE EXCHANGER 
BAFFLES: - 


Single Segmantal with cut Cii of Shell IDD 

Baffle Shell Diametral Clearenee CmnO 

Baffle pitch at Inlet End CnatD 

Baffle pitch at Outlet End CmnO 

Baf f leHole-TubeOD Diametral Clearenee CmnO 
Net window flow area Csq_HD 

Minimum cross flow Area Csq_nO 

Number of Pairs of Sealing Strips C^ 

Cross flow Reynolds number C;iO 


= 25 
= 3.81 
= 375 
= 375 

= 7.9375E-1 
= 2. 7337E-2 
= 2-8575E-2 
= 2 

= 2.0575E+4 


PRESSURE DROP 

Correction factors, Rl 

Rb 

Rs 

Cross Flow Pressure Drop 


CN/sq_TiD 

CN/^sq_nO 


4.4S26E-1 
1 . 9444E-1 
2. 401SE-1 
1. 1635E-I-4 
1 . 2810E-*-4 


Window Pressure Exrop 


End zones Pressure EJrop 
Total Shell side Pressure drop 
HEAT TRANSFER COEFFICINT 

Ideal C Cross flowU Heat Trans. Coeff 
Correction factors* Jc 

J1 
Jb 
Js 

Cummulative correction factor 
Shell Side Heat Transfer Coefficient 

VIBRATION CHECK REPORT 
For Tubo m&tAl d&nsity = 7S4^ kg/ m 

Modulus of elasticity = a.OEll N^" 

Fatigue stress = 3. OES N/ m 

For Span Length « Central Baffle Pitch 
Fundamental tube natural frequency = 4. 6SS539SSS9E+OS Hz 

Vorte>c shedding Frequency = S . SBl 7SSSS4i E+OO Hz 

Critical cross flow velocity = 3. 7394935041 E+Ol Cm/s!> 

Maximum cross flow velocity = S. 733a494436E-01 Cm^sD 

Baffle damage number = 4. 355533S994E— 03 

Col 1 i s i on damag e number = S . 41 01 1 91 51 OE— 04 

Maximum damage number = 5. S7937391 43E-01 

Safe in Vortex shedding 

fvs/fn = 1.761 131 9361 E-Oa < 0.5 

Safe in fluid elastic whirling 

vcrit = a. 7394965541 E+01 > vmax = 6. 733a494436E-01 
Safe No Baffle Damage 

Nbd = 4. 355533S994E-03 < Max allow number = 5. 87937391 43E-Oi 

Safe ««« NO Collision Damage 


S3 

CN^^'s g g p = 3. 41S0E+3 
CN^sq_mD = 3. SS64E+4 

. CWy^'sq_m. K> = 8. 4030E+3 
= 1.0340 
= 6.S96E-1 
= 9.asaE-i 
= 9.siaE-i 
= 6. 133E-1 
CW/sq_m. KD = 5. 14BE+3 


S3 


Ned = S. 4,101191610E-04 < Max allow number = B. 87937291 43E-01 


For- Span Length ■ 2 # C Baffle Pitch! 

Fundamental tube natural frequency = 1 . 1 71 3S23972E-t-02 Hz 

Vortex shedding Frequency = S. 2B17SSS941E-*-00 Hz 

Critical cross flow velocity = 6. S4S74S3SB2E-<-00 Cnv's! 

Maximum cross flow velocity = S. 7332494426E-01 CmXs! 

Baf f 1 e damage number = 1 - 74221 3BB97E— 02 

Collision damage number = 1 . 34BS1 90S42E-02 

Maxi mum damage number = B. S7927291 43E-01 

Safe in Vortex shedding 

fvs/fn = 7. 0444S77042E-02 < O. B 

Safe ««« in fluid elastic whirling 

vcrit = e. S4S74e2SB2E-»-00 > vmax = S. 733249442SE-01 
Safe No Baffle Damage 

Nbd = 1 . 74221 3BB97E -02 < Max allow number = B. S792729143E-01 

Safe NO Collision Damage 

Ned = 1 . 34B61 90642E-02 < Max allow number = B. S792729143E-01 


For Span Length ■ Max. urisupported length 
C First window from ends! 


Fundamental tube natural frequency = 7. 49SS473422E-»-01 Hz 


Vortex shedding Frequency 
Critical cross flow velocity 
Maximum cross flow velocity 
Baffle damage number 
Collision damage number 
Maximum damage number 


= S. 2B17SSS941E-t-00 Hz 
= 4. 3S3i97SS6SE-i-00 Cm/s! 

= S. 733249442SE-Oi Ciiv''s! 
= 2 . 72220SS871 E-02 

= 3. 2SB2027933E-02 

= B. 27927291 43E-01 


Safe in Vortex shedding 

fvs/'fn = 1 . 100701 203SE-01 < O. B 

Safe in fluid elastic whirling 

vcrit = 4 . 3S31976S66E-1-00 > vmax = S. 7332494426E-01 

Safe ««« No Baffle Damage 

Nbd = 2 - 7222086871 E-02 < Max allow number = B. 87927291 43E-01 

Safe NO Collision Damage 

Ned = 3 . 28S2027933E-02 < Max allow number = B. 87927291 43E-01 



